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ABSTRACT 

We present the results of a new Hi, optical, and Ha interferometric study of the 
nearby spiral galaxy NGC 157. Our combined C- and D-array observations with the 
VLA show a large-scale, ring-like structure in the neutral hydrogen underlying the 
optical disk, together with an extended, low surface density component going out to 
nearly twice the Holmberg radius. Beginning just inside the edge of the star-forming 
disk, the line of nodes in the gas disk commences a 60° warp, while at the same time, 
the rotation velocity drops by almost half its peak value of 200 km s~^, before leveling 
off again in the outer parts. While a flat rotation curve in NGC 157 cannot be ruled 
out, supportive evidence for an abrupt decline comes from the ionised gas kinematics, 
the optical surface photometry, and the global H i profile. A standard 'maximum- 
disk' mass model predicts comparable amounts of dark and luminous matter within 
NGC 157. Alternatively, a model employing a disk truncated at 2 disk scale lengths 
could equally well account for the unusual form of the rotation curve in NGC 157. 
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1 INTRODUCTION 

Observations of the neutral hydrogen in 'normal' spiral 
galaxies have frequently revealed new and surprising as- 
pects of their structure. Structures in the gas such as warps, 
tails, and superbubbles point to a more turbulent history 
than that indicated by the stellar distribution. The nearby 
SAB(rs)bc galaxy NGC 157 is a striking example of this. 
Blue light photographs (Lynds 1974) show a dusty, floc- 
culent disk of high surface brightness, which is reasonably 
symmetric and apparently undisturbed. In contrast, the 
Hi aperture synthesis observations presented here show a 
kinematic warp and an unusual rotation curve, which may 
point to a strongly warped disk or a peculiar dark matter 
distribution. 

Some hints of NGC 157's distinctive qualities are con- 
tained in earlier optical studies. Blackman (1979) presented 
photographic UBVR luminosity profiles, which pointed to 
the presence of a second, outer exponential disk, with a 
longer scale length than the inner disk. He also conjectured 
(from modeling of the spiral arm morphology) that the ro- 
tation curve might be declining at large radii. The disk kine- 
matics of NGC 157 were first investigated by Burbidge, Bur- 
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bidge & Prendergast (1961), who derived a slowly-rising ro- 
tation curve and a total mass ~ 6 x 10^" M0 on the basis of 
two long-slit spectra close to the galaxy's major axis. Zasov 
& Kyazumov (1981, hereafter ZK) obtained long-slit spec- 
tra of NGC 157 at a number of position angles, and found 
good agreement with the Burbidge et al. rotation curve in 
the southwest part of the disk, but a significant drop in 
rotational velocity (up to ~ 100 km s~^) beyond a radial 
distance of 55 arcsec to the northeast. Further evidence of 
peculiar kinematics comes from the global Hi profiles of 
NGC 157 obtained with the Jodrell Bank 76-m radio tele- 
scope by Staveley-Smith & Davies (1987; hereafter SSD), 
and with the Parkes 64-m telescope by Mathewson, Ford & 
Buchhorn 1992; hereafter MFB). Both spectra show a two- 
horned profile with a strange pair of shoulder- like features. 
All this is made the more curious by the fact that the galaxy 
is quite isolated. No close companions of comparable bright- 
ness are seen in the POSS images, and the nearest galaxy in 
the NED database within 5{cz) = ±500 km s"^ is UGCA 14, 
at a projected distance of 3?5 (~ 1.3 Mpc for D = 20.9 Mpc; 
TuUy 1988). In fact, there are only five catalogued galaxies 
within 5° and 5(cz) — ±1500 km s~^. 

There are also hints from the literature of strong inter- 
nal forces acting on the galaxy. At the position of the pecu- 
liar motions seen by ZK, blue light photographs show no ob- 
vious signs of a perturbation or companion object. However, 



© 1997 RAS 



2 S. Ryder, A. Zasov, V. Mclntyre, W. Walsh and O. Sil'chenko 



Table 1. VLA observing parameters for NGC 157. 
Parameter Value 



a (J2000 pointing centre) 
5 (J2000 pointing center) 
Distance adopted 
Natural-weighted beam FWHP 
Uniform-weighted beam FWHP 
Channel map rms noise (natural) 
Channel map rms noise (uniform) 
Continuum map rms noise (uniform) 
IF 1 central frequency 
IF 2 central frequency 
Channel increment (unsmoothed) 



00" 34"" 46! 5 
-08° 23' 48" 
20.9 Mpc 
41" X 27" (4.2 X 2.8 kpc) 
18" X 12" (1.8 X 1.2 kpc) 
0.6 mjy beam~^ 
1.3 mJy beam~^ 
0.3 mJy beam~^ 
1412.9 MHz 
1411.6 MHz 
5.2 km s-i 



Hq images (ZK; Rozas, Beckman & Knapen 1996) reveal a 
'ring' of H ll regions centered ~ 1 arcmin northeast of the 
nucleus of NGC 157 and some 40" (~ 4 kpc) in diameter. A 
ring-shaped feature like this could be caused by propagating 
star formation on the periphery of an expanding H I 'super- 
bubble'. If so, then (relative to its parent galaxy) it would 
be one of the largest and most energetic examples found to 
date. 

In order to unravel some of these mysteries, and to gain 
a better understanding of the morphology and evolution of 
NGC 157, we have made the first Hi map of NGC 157; ob- 
tained new wide-field B, I, and Ha imaging; and mapped the 
kinematics in detail of the inner disk using Ha Fabry-Perot 
interferometry. Details of these observations, and their re- 
ductions, are given in Section |^. An analysis of the galaxy's 
kinematical and morphological parameters is presented in 
Section Our interpretation of these new results is dis- 
cussed in Section ^, and a summary of our conclusions fol- 
lows in Section |^. 



2 OBSERVATIONS 
2.1 Hi Mapping 

NGC 157 was observed with the two lowest-resolution con- 
figurations of the Very Large Array (VLA) of the National 
Radio Astronomy Observatory (NRAO). A total of 3.7 hours 
on-source integration time was obtained with the C-array on 

1994 October 27, and a further 2.8 hours with the D-array on 

1995 May 26. The primary fiux density (and bandpass) cali- 
brator was 3C48, assumed to have a fiux density of 16.12 Jy 
at 1.413 GHz. The nearby radio source 0022-1-002 (J2000) 
was used as a secondary calibrator to monitor variations in 
the gain and phase. To be sure of covering the full projected 
velocity width of NGC 157 (324 km s'^ SSD) at modest 
resolution (10.4 km s~^ after on-line Hanning smoothing), 
the two Intermediate Frequencies (IFs) were configured to 
64 channels each with 1.56 MHz bandwidth. By tuning the 
IFs to two closely-spaced frequencies (Table |l]) spanning that 
of the redshifted H i line, it was possible to accommodate the 
full profile width, with a few channels of continuum at either 
end of the band, and with 8 channels of overlap between the 
IFs. 

The C-array and D-array UD-datasets were processed 
separately using implementations of the NRAO aips soft- 
ware at the NRAO Array Operations Centre and at the Aus- 
tralia Telescope National Facility. The visibility data were 



calibrated using standard VLA procedures, and bad data re- 
moved interactively. After combining the two datasets in the 
uu-plane, the imaging task MX was used to construct a cube 
of dirty channel maps for each of the two IFs, as well as cubes 
of the corresponding dirty beam shapes. At this stage, the 
two cubes were 'sliced' into their constituent channel maps, 
and then a 'supercube' was constructed by stacking all avail- 
able maps in order of increasing frequency. After dropping 
some channel maps at either end of the IF bandpass with 
low signal-to- noise ratio, we were left with 112 usable chan- 
nel maps in the supercube spanning 581 km s"^. 

A total of 20 channels to either side of the main pro- 
file containing only continuum emission were averaged to- 
gether to form a master dirty continuum image, and then 
this image was subtracted from every map in the supercube. 
Although continuum subtraction in the uv-plane is always 
preferable over subtraction in the image domain (Killeen 
1993), the need to combine two IFs made this impractical. In 
any case, since Au/i/q ~ 0.002, and the continuum emission 
from NGC 157 itself is not particularly strong, this method 
of continuum subtraction worked extremely well. 

Because they highlight different structural features of 
NGC 157, two distinct supercubes were constructed using 
different weighting schemes in the imaging process. Uniform 
weighting yields the best possible spatial resolution, while 
natural weighting offers improved sensitivity to extended, 
low surface brightness emission, both of which turn out to 
be important in understanding the nature of NGC 157. In 
each case, the pixel scale was chosen to allow at least 3 pix- 
els across the synthesized beam minor axis (Table |^), and 
the map size to match the VLA primary beam at 1.4 GHz 
(~ 30 arcmin). By comparing the total CLEANed flux den- 
sity with the improvement in rms noise obtained, as a func- 
tion of the number of CLEAN iterations performed using 
the APCLN task, we elected to perform 200 CLEAN itera- 
tions on each map of the natural-weighted supercube, while 
the uniform-weighted data was found to require no further 
CLEANing. 

Following primary-beam corrections, and conversion 
from frequency to optical heliocentric velocity, moment 
maps (H I column density, velocity, and velocity width) were 
constructed from each of the two supercubes, using masks 
derived from a spatially and spectrally smoothed version of 
the supercube, and by setting a fiux cutoff that just enabled 
a clean separation of the H I emission from background noise. 
Beam full-widths at half power (FWHP) and rms noise lev- 
els in the final CLEANed maps are summarised in Table ^. 



2.2 Optical Imaging 

NGC 157 was observed under photometric conditions with 
the Kitt Peak National Observatory 0.9-m telescope on 
1995 July 3 UT. A Tektronix 2K x 2K CCD was used at 
//7.5 with the doublet corrector, making available a field of 
23.2 arcmin x 23.2 arcmin at 0.68 arcsec pixel"^. Two 500 s 
exposures in B, and two 300 s exposures in I were obtained, 
together with short bracketing exposures of Landolt's (1983) 
Field 109. On the previous night, NGC 157 was observed 
once for 1000 s in the KP1564 (6618 A central wavelength, 
74 A FWHM) Ha filter, and also for 1000 s in the KP808 
(6411/88) red continuum filter. Images of Cygnus OB2 No. 9 
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(Massey et al. 1988) with the KP1564 fiher served as a flux 
calibration for the Hq data. 

Following removal of the bias level and structure, at- 
tempts were made to flatfield the data using dome flatfields, 
twilight sky flatfields, or both, but residual structure at the 
5 per cent level was still apparent near the edges of the 
frames. Finally, low-order polynomial fits to these border re- 
gions were computed and used to normalise both the object 
and calibration frames. Aperture photometry of the broad- 
band standards over a large range in airmass allowed us 
to compute photometric zero points, extinction, and colour- 
correction terms. Each pair of B and / images of NGC 157 
were averaged together, after scaling to a common modal 
value in the central regions, and registration using field stars. 
The red continuum image of NGC 157 was convolved to the 
same stellar FWHM as the Hq image, registered using field 
stars, multiplied by a factor of 1.2 (the ratio of the contin- 
uum to line filter full widths) , and then subtracted from the 
KP1564 image, leaving just the pure Ha emission. 



2.3 Ha Fabry- Perot Interferometry 

The two-dimensional velocity field for the ionized gas in 
NGC 157 was obtained at the 6-m telescope at the Special 
Astrophysical Observatory of the Russian Academy of Sci- 
ences on 1995 October 24, when the seeing averaged 2.5 arc- 
sec. A scanning Fabry-Perot interferometer was installed in 
the pupil plane of a focal reducer which was attached to the 
//4 prime focus of the telescope (yielding an //2.4 output 
beam). The detector used was a 512 x 512 pixel intensifled 
photon-counting system. The interferometer was operated in 
the 501st order at 6562.8 A, and a narrow-band (10 A) fil- 
ter at the appropriate Ha redshift served to block adjacent 
orders. A total of 32 images (each 4 arcmin x 4 arcmin) 
were obtained with a velocity increment of 18.84 km s~^, 
with 3 minutes of integration at each velocity step. A neon 
emission line at 6598.95 A was used for wavelength calibra- 
tion. The raw observational data were rebinned to cubes 
of 32 X 256 X 256 pixels, with an image scale of 0.92 arc- 
sec pixel"^ and a spectral resolution of 2 — 2.5 channels 
(40 — 50 km s~^). Following the steps outlined in Le Coarer 
et al. (1992) and Laval et al. (1987), standard reduction pro- 
cedures (corrections for phase shifting, subtraction of the 
night-sky emission-line spectrum, construction of the veloc- 
ity field, etc.) were performed using the ADHOC Fabry- 
Perot software package developed at the Marseille Observa- 
tory (Boulesteix 1993). 



3 RESULTS AND ANALYSIS 
3.1 H I distribution 

The resulting channel maps from the natural-weighted su- 
percube are shown in Figure 0. Already, there are indica- 
tions from this sequence of maps of some peculiarities in 
the Hi structure and dynamics; namely the existence of an 
extended gas disk beyond the optical radius (r25 — 2.08 ar- 
cmin; de Vaucouleurs et al. 1991, hereafter RC3) of this 
galaxy, and from the bending of the isovelocity features, an 
associated warping of the gas disk. 

Before considering this disk in more detail, it is prudent 



to examine the extent to which we may be losing flux and 
large-scale structure due to missing baselines. To do that, we 
have computed the global H I profile for NGC 157 by using 
a moving window to optimise the emission measured within 
each channel map of the natural-weighted supercube, and 
plotted the fluxes versus the respective map velocity in Fig- 
ure^ Our profile matches well the general form of those pub- 
lished by SSD and by MFB, particularly the 'shoulders' at 
~ 1540 km s^^ and ~ 1800 km s~^, although our peak mea- 
sured fluxes are slightly higher. A comparison of our profile 
moments (integrated flux, W20, etc.) with those from SSD 
and MFB appears in Table ^ We take the good agreement 
shown there as evidence that our VLA observations are not 
resolving out a signiflcant low column density component, 
at least within the 12 — 15 arcmin single dish beamwidths. 

The shape of the spectrum can be shown to be the re- 
sult of the superposition of two spatially distinct compo- 
nents in the gas distribution. The dotted and dashed lines 
in Figure ^ are spectra that were obtained by separately 
summing the regions inside and outside, respectively, of the 
fiB = 25 mag arcsec"'^ isophote. A similar result is ob- 
tained if the summation is done using a mask chos en t o 
separate the 'warped' component of the Hi (Section 3.2). 
The flux ratio of the two components is (inner/outer) = 
(37±2 Jy km s"^)/(41±2 Jy km s"^) = 0.9. Their systemic 
velocities differ slightly, being 1674 ±3 and 1685 ±3 km s"'^ 
respectively. 

A flux integral of (78 ± 4) Jy km s~^ corresponds to a 
total mass of Hi of (7.9 ± 0.4) x 10^ Mq for the distance 
of 20.9 Mpc (TuUy 1988) adopted here. Taking the abso- 
lute blue magnitude for NGC 157 to be Mb = -20.8 (Sec- 
tion ^.5[ ), we therefore find Mm/ Lb ~ 0.23, quite typical of 
an Sbc-type galaxy (Roberts & Haynes 1994). 

The Hi column density distribution (zeroth moment) 
is shown as contours superimposed upon a section of the 
KPNO 0.9-m _B-band image in Figure ^, and as a grey- 
scale in Figure ^. The impression of an extended Hi disk 
in NGC 157 gained from the channel maps is confirmed 
by these images, with the outermost contour, represent- 
ing a projected H l column density of 0.5 x 10^" cm~'^ 
(0.4 Mq pc~^) spanning some 11 arcmin, or over 2.5 times 
the optical diameter. The flux cutoff imposed in defining the 
edge of the Hi disk is roughly at 3x the noise level of the 
moment fitting; however, in blanking off the low signal-to- 
noise regions, we were guided also by the smooth continuity 
of the velocity field in the outer disk (Figure ^). NGC 157 
possesses a central Hi 'hole' of the sort commonly seen in 
spiral galaxies, and the projected gas column density reaches 
a peak in a ring underlying the spiral arms. An arm or tail- 
like feature appears to be breaking free from the northern 
edge of the H l disk in Figure ^. 

The apparent offset between the major axis position 
angles of the optical and gas disks in Figure ^ seems to be 
the result of a change in the orientation of the disk orbits 
between the inner and outer disk. The isovelocity contours in 
Figure ^ show the kinematic line-of-nodes precesses through 
~ 60° between the edge of the optical disk and the limits of 
the outer gas disk. At this resolution, the inner disk displays 
the characteristic pattern for an inclined disk in differential 
rotation, albeit with an abrupt turnover in the p rojec ted 
velocity as a function of radius (discussed in Section |4.3| ) . In 
the outer disk, the isovelocity contours meander in a fashion 
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Table 2. Global Hi profile parameters for NGC 157. 

Parameter Staveley-Smith & Davies 1987 Mathcwson ct al. 1992 This Study 

Flux Integral (Jy km s"!) 62.6 ± 5.5 72.8 ± 7.1 78.0 ± 4.0 

Vhci" (km s-i) 1669 ± 7 1655 ± 7 1682 ± 12 

W20 (km s-i) 324 ± 14 336 ± 14 326 ± 10 

"Mean of velocities where the H I profile falls to 20 per cent of the peak flux (Staveley-Smith 
& Davies 1987; this study), or 50 per cent of the peak flux (Mathewson et al. 1992) 



typical of galaxies with a warp. There are also kinks in these 
contours consistent with the action of a density wave as 
mentioned earlier. 

A higher resolution map made with uniform-weighting 
(Figure ph resolves out nearly all of the extended outer H I 
disk, but does highlight some interesting structure within 
the inner disk. The inner gas 'ring' is now revealed to be 
composed of a honeycomb-like lattice of rings and cells. Al- 
though initially chaotic in appearance, the peaks in the gas 
surface density making up the ridge-lines of these 'cells' are 
shown to be intimately associated with the spiral arms and 
other sites of massive star formation, when contours of the 
Ha surface brightness are superimposed. The Ha ring to the 
northeast of the galaxy discussed by ZK is seen to have a 
distinct H I counterpart, as do many of the other major star- 
forming complexes in NGC 157. At the highest resolution 
offered by the Hi data, the individual cells have diameters 
between 2 and 4 kpc, and deprojected Hi column densities 
peaking at ~ 13 M© pc~^ on the cell 'walls', dropping to as 
low as 1.5 Mq pc~^ at the centres. The correlation of star 
formation activity with local gas surface density in NGC 157 
is somewhat better than that typically observed in spiral 
galaxies (e.g., Ryder et al. 1995; Rownd et al. 1994), and is 
more like that observed in the shells of the dwarf irregular 
galaxy Holmberg II (Puche et al. 1992), although there is 
no suggestion that the gas 'cells' in NGC 157 are expanding 
coherently. 

The azimuthally-averaged column density of H I as a 
function of radius in the naturally-weighted zeroth moment 
map is plotted in Figure ^. Since the intention here is to 
indicate the column density of gas in the disk of the galaxy 
that would be seen from 'face-on', we have adopted a fixed 
incl ination (i = 45° based on the results of Sections 3.2 and 
as well as Grosb0l (1985)) and position angle 
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220°) 

and then calculated the deprojected column densities ac- 
cordingly. The arrow in Figure M indicates the radius corre- 
sponding to the fiB ~ 25 mag arcsec"^ isophote (^25), and 
just on half of the total Hi content of NGC 157 lies outside 
this radius. The point at which the Hi surface density falls 
below 1.0 M0 pc~^ equates to 2.1r25, an H l-to-optical ex- 
tent bettered by only 5 of the 23 comparable galaxies studied 
by Broeils & van Woerden (1994; hereafter BvW). 

3.2 H I kinematics 

We determined a rotation curve from the velocity field by 
using an AlPS implementation of the ROTCUR algorithm 
(Begeman 1989) to fit rings of variable inclination and posi- 
tion angle under the assumption of uniform circular mo- 
tion in each ring. To begin with, we analysed only the 
regions within the optical extent of NGC 157 using our 
highest resolution (uniform-weighted) velocity field, sam- 



pled at approximately one synthesized beamwidth (15 arc- 
sec) intervals. By adopting starting values for the kinemat- 
ical parameters (dynamical center position, systemic ve- 
locity, and disk tilt and position angle) d raw n from RC3 
and our own surface photometry (Section 3.5), and allow- 
ing all of these parameters to vary freely, we arrived at the 
following values for the dynamical center: Vhd = 1676 ± 
9kms-\ adyn(J2000) = 00*^34"" 46':7(±4"), 5dyn(J2000) = 
-08° 23' 49" (±3"). These results are in close agreement 
with the position of the optical centre (shown as the point- 
ing centre in Table ^ and the systemic velocities from the 
midpoint of the global Hi profile (Table |^. 

The analysis was then repeated with the dynamical 
center held fixed at these values, so that the combinations 
of rotational velocity, inclination, and position angle (kine- 
matic line-of-nodes) which minimised the systematic veloc- 
ity residuals within each ring could be studied. Fixing the 
dynamical centre also allows us to compare the results from 
fitting the entire disk with those that come from fitting only 
the 'receding' side (i.e., mostly the southern half) or the 'ap- 
proaching' (northern) side, which then highlights any asym- 
metries in the disk kinematics. This same procedure can 
be applied to the full H I disk kinematics as revealed by the 
natural-weighted data, although at a somewhat coarser reso- 
lution (35 arcsec). The dynamical centre in this case is found 
to lie within 2 arcsec and 1 km s~^ of that found earlier, so 
has been fixed at identical values. 

The end result of both of these analyses is shown in 
Figure ^, in which points within 120 arcsec of the centre are 
drawn from the uniform-weighted fit, and those further out 
are from the natural-weighted data. Clearly, the continuity 
between these two is extremely good, and the major differ- 
ences are in fact between the two halves of the disk. The 
principal feature of the rotation curve is a strong decline 
that begins just inside r2s. The exact shape of the decline 
is somewhat uncertain because the fits to the northern (re- 
ceding) half failed to converge in the transition region, and 
the two halves of the disk give rather different results for 
r = 200 — 300 arcsec. However, the fits to the full disk and 
southern half do agree in the transition region, and all three 
fits agree for r > 350 arcsec. At this point, the circular veloc- 
ity seems to have settled to a constant value at ~ 120 km s~^ 
(55% of the maximum rotational velocity), and the inclina- 
tion is only slightly higher than that attained in the inner 
disk. 

Although the steady progression in position angle is 
characteristic of a warp, the rotational velocity drops by al- 
most half in less than 1 arcmin, while the inclination varies 
by no more than 10° in this region. The drop in velocity is 
even faster than Keplerian for 120 < r < 190 arcsec, though 
such be haviou r is permissible for a disk with a sharp edge 
(Section EZi). 



© 1997 RAS, MNRAS 000, §-26 



The Peculiar Rotation Curve of NGC 157 5 



3.3 Ionised Gas Kinematics 

Figure |^ shows a comparison of the Ha Fabry-Perot veloc- 
ity field with the highest resolution velocity field available 
from the uniform- weighted Hi data. Owing to the differ- 
ing resolutions, and the complexity of the Ha velocity field, 
only selected velocity ranges and isovelocity contours for the 
Hi data are shown. There appears to be a discrepancy be- 
tween the velocity scales of perhaps 10 — 20 km s~^, but 
this is within the spectral and spatial resolutions of the two 
datasets. Otherwise, the agreement between the neutral and 
ionised gas kinematics in tracing structure within the inner 
velocity field is remarkably good. 

A rotation curve was derived from the two-dimensional 
velocity field using the ROTCUR algorithm described in 
Section B.2. Preliminary fits indicated a rotation centre 
again coincident with the optical nucleus, but a systemic 
velocity of 1671 ± 7 km s~^, consistent with the velocity 
offset just mentioned and still within the range of Hi sys- 
temic velocities in Table ^. Keeping these parameters fixed, 
the best-fitting combination of rotation velocity, inclination 
and position angle have been evaluated, and plotted as solid 
lines in Figure M. Unusual structure in the nuclear veloc- 
ity field, possibly associated with a bar (Section 4.2), makes 
fitting the rotation curve there difficult, but in the range 
30 < r < 70 arcsec, consistency between the Ha, Hi, and 
photometric orientation parameters is excellent. 

The Ha rotation curve rises relatively slowly, as a higher 
inclination is initially favoured, but then accelerates as the 
inclination drops toward the 50° angle of the inner disk. It 
then levels off at the same peak velocity (~ 200 km s"'^) at- 
tained by the H I. A bulge component may be responsible for 
the steep central velocity gradient. It is worth noting, how- 
ever, that the inner maximum of the rotation curve cannot 
be accounted for by a combination of a King-profile bulge 
(even if it is sharply truncated) and the disk, but can be fit 
by the combination of a steep, inner expon ential disk and a 
flatter, outer disk (Blackman 1979; Section [iX^. 



3.4 1.4 GHz continuum 

Since only one of the IFs contained a significant number of 
line-free channels, the AlPS task MX was used for both imag- 
ing and CLEANing of the 1.4 GHz continuum (in prefer ence 
to the continuum map yielded as a by-product in Section ^.l[ 
which was formed by averaging in the image plane of non- 
contiguous channels, and which could only be CLEANed 
over a limited area). A total of 17 channels was extracted 
from IF 2 of the -uw-dataset, Fourier transformed to the im- 
age domain using uniform-weighting, and then subjected to 
1000 iterations of the CLEAN algorithm using a similarly- 
extracted dirty beam. Contours of the 1.4 GHz continuum 
emission have been superimposed on top of the Ha image of 
NGC 157 in Figure |. 



As was the case with the Hi (Section 3.1), there is a 



close, though not strict correspondence of the radio contin- 
uum fiux density with the Ha surface brightness. A contin- 
uous ridge of emission underlies both the major spiral arms, 
passing through the nucleus and reaching a peak midway 
along the southern arm, but with an abrupt break shortly 
thereafter. Diffuse emission is seen to fill the entire opti- 
cal disk, but does not extend much beyond this. All of this 



is consistent with the origin of the continuum emission be- 
ing mostly from relativistic electrons diffusing along galactic 
magnetic field lines, but with some contribution from the 
thermal component from Hll regions. Interestingly, the ra- 
dio continuum contours conform also to the shape of the 
northern ring seen in Hi and in Ha. 

An unresolved radio continuum source of some 146 mjy 
lies just over 5 arcmin southeast of the nucleus of NGC 157. 
It has no counterpart in existing radio or optical catalogues. 
The source lies well outside the optical disk, but is just in- 
side the detectable limit of the H I disk. Because of the high 
galactic latitude of NGC 157 {b — —71°), this object is al- 
most certainly at a greater redshift than NGC 157. This 
opens up the tantalising possibility of measuring H I absorp- 
tion, and thus the gas spin temperature, in a part of the gas 
disk where such measurements are rare. 

We have searched for signs of such absorption in both 
the full uu-dataset, and in the continuum-subtracted cube, 
but to no avail. Such a null result is perhaps to be expected, 
given the conclusion of Dickey, Brinks & Puche (1992) that 
even the longest baselines of the VLA C- and D-arrays (3 km 
maximum) are inadequate for resolving out the H I emission 
fiuctuations across the continuum source. FoUowup observa- 
tions with the VLA A- or B-array may yet however succeed 
in placing limits on the gas temperature in the very outer 
parts of the gas disk in NGC 157. 

3.5 Surface photometry 

The GASP software package was used to compute radial sur- 
face brightness profiles and disk orientation parameters for 
NGC 157 in a manner similar to that described by Ryder 
& Dopita (1994) for a large sample of southern spiral galax- 
ies. The B, I, and continuum-subtracted Ha images were 
first block-averaged to 1.36 arcsec pixel"^, after which the 
cores and extensive halos of the two bright foreground stars 
visible in Figure ^ were masked out. The sky background 
was determined from the modal peak of the histogram of all 
data values within 20 pixels of the image border, and found 
to be 22.01 ± 0.02 and 19.28 ± 0.01 mag arcsec"^ in B and 
I respectively. 

For a given semi-major axis, GASP iteratively attempts 
to fit ellipses to the galaxy image, varying the ellipticity and 
position angle (while holding the ellipse centre fixed at the 
position of the nucleus) so as to minimise residuals about the 
mode of the pixel values around the ellipse perimeter. The 
semi-major axis is increased by 15 per cent each time (the in- 
creased number of pixels sampled at each radius partly com- 
pensating for the reduced signal-to-noise in the outermost 
pixels), until a level equal to the measured sky background 
is reached. 

The modal surface brightnesses from separate B and 
/ analyses, corrected for atmospheric (but not Galactic or 
internal) extinction, are plotted in Figure ^ Also shown 
on this plot is the mean surface brightness in Ha (with 
-1-24.0 mag arcsec"^ being equivalent to an Ha flux density 
of 10~^^ ergs cm~^ s~^ arcsec"^) computed at each radius 
using the ellipse parameters from fitting to the /-band im- 
age. The variation in ellipse inclination and position angle 
as a function of radius in the I band (less affected by dust 
than B) has also been plotted in Figure |^ for comparison 
with those determined from the gas kinematics. 
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Blackman (1979) reported the presence of an extra lin- 
ear outer component, in addition to the inner bulge and 
exponential disk. His reduced luminosity profiles are not di- 
rectly comparable to our radial surface brightness profiles, 
but Figure [lo| does confirm the existence of a second ex- 
ponential disk in the range (90 ^ r ^ 200 arcsec), with a 
longer scale length than the inner exponential disk, which 
covers (35 ^ r 90 arcsec). This outer exponential disk also 
coincides with the transition region between the inner and 
outer H l disks, as well as the region of most rapid turnover 
in the rotation curve. 

The integrated intensity out to the sky background 
yields apparent magnitudes in B and / of 11.29 ± 0.03 and 
9.47 ± 0.02 respectively. After corrections for Galactic ex- 
tinction (As — 0.12; Burstein & Heiles 1984) and for ex- 
tinction internal to NGC 157 at an inclination of (45 ± 5)° 
(TuUy & Fouque 1985), we arrive at corrected total mag- 
nitudes (in the RC3 convention) of = 10.79 ± 0.06 and 
It = 9.28 ±0.05. Adopting the same distance (20.9 Mpc) as 
used for the total H I mass determination yields an absolute 
magnitude Mb ~ —20.82 ±0.06, and thus a total blue lumi- 
nosity of ~ 3.3 X 10^" Lq. Our B^ value is consistent with 
equivalent tabulations in RC3 and TuUy (1988), although 
the individual magnitudes and extinction corrections do dif- 
fer. MFB found some 0.3 mag fainter than us, though 
the integrated I magnitudes are in much better agreement, 
so much of the difference lies in the (still contentious) cor- 
rections for extinction. 

The continuum-subtracted Ha data was similarly inte- 
grated out to D25, resulting in a total flux in the Ha line 
of (2.1 ± 0.2) X 10"^^ ergs cm'^ s"^ (~ 2.7 x 10® Lq for 
the adopted distance). This includes a correction for atmo- 
spheric extinction, as well as the 1.1 mag of internal ex- 
tinction applied by Kennicutt (1983) in his study of the 
star formation rate (SFR) in normal disk galaxies. Using 
the same conversion formulae from Ha luminosity to SFR 
as Kennicutt, we estimate a total SFR for massive stars 
(M ^ 10 M0) of 1.6 Mq yr~^, and over all stellar masses 
(0.1 Mq < M < 100 Mq) to be ~ 10 Mq yr~^ Nor- 
malising yields the following (distance-independent) quan- 
tities: (SFR/disk area) = 30 Mq pc'^ Gyr-^; (SFR/Lb) = 
0.3 Mq Gyr-i L'^; and (SFR/Mm) = 1.3 Gyr'^ Con- 
sequently, although the absolute present-day SFR is not 
unusual for its type (Sbc), relative to its size and mass, 
NGC 157 is nearly twice as active in forming stars as galax- 
ies of similar type and colour (Ryder 1993). 

This enhancement of Ha emission is conflrmed by the 
far-infrared to Ha luminosity ratio, Lfiu/Lhc- The value 
of logLpiR, obtained from IRAS data after reducing to 
the adopted distance is 10.34 (Young et al. 1989); hence 
logLpiR/LHa ~ 1.9, whereas for bright galaxies of a similar 
colour, typical ratios are in the range 2.1 — 2.5 (Zasov 1995). 
Note that log Lpm/ Lb = -0.18, observed in NGC 157, is 
quite normal for late-type spiral galaxies (Young et al. 1989; 
Zasov 1995), in contrast to the normalised Ha star forma- 
tion rates. Using the relationship between star formation 
rate and far-infrared luminosity (SFR ~ 2.5 x IO^^I/fir/Lq; 
Zasov 1995), we obtain SFR ~ 5A/q yr~\ No Ha emis- 
sion was detected in the extended Hi disk, to a Scr limit of 
3 x 10~^^ ergs cm~^ s"'^ arcsec"^, implying that the mas- 
sive star formation rate (M J; 10 Mq) there cannot exceed 
- 2 Mq pc"^ Gyr-\ 



4 DISCUSSION 

4.1 A superbubble in NGC 157? 

Given that one of the primary motivations for this study 
was the possible existence of a major H I superbubble in the 
disk of NGC 157, how convincing is the new evidence for or 
against such a feature? Normally, the strongest evidence for 
such superbubbles comes from the localised distortions they 
introduce to the global velocity field (Ryder et al. 1995) , and 
yet the isovelocity contours in Figure ^ (as well as the higher 
resolution images) exhibit a remarkable degree of symmetry 
about both the major and minor axes. Examination of the 
position-velocity diagrams in this region also indicates no 
unusual profile splitting or deviations attributable to the 
action of an expanding superbubble, although the fact that 
the projected rotation velocity is changing rapidly in this 
particular region complicates such an analysis. 

How then are we to explain the abrupt change in ve- 
locity across the northeast ring seen by ZK in Ha? We have 
made cuts through the uniform- weighted Hi cube at simi- 
lar position angles (e.g., d = 10°; Figure |ll|), and find that 
beyond 45 arcsec along the northeastern (approaching) axis, 
there is indeed a drop of almost 40 km s~^ in velocity which, 
when deprojected by the same inclination (i = 31°) used 
by ZK and Burbidge et al. (1961), results in a total drop 
of ~ 80 km s~^, as reported by ZK. The Hi cut however 
shows a similar drop on the opposite axis beyond 60 arcsec, 
by which point ZK were unable to detect any more optical 
emission. Thus, the velocity discontinuities reported by ZK 
are real, but are caused by something other than the action 
of an H I superbubble. 



4.2 The bar in NGC 157 

Surface photometry in the near infrared {JHK) bands has 
recently been published by Elmegreen et al. (1996a, b), who 
make the case for a bar component, with a flat luminosity 
profile going out to one-quarter of r25. Sempere & Rozas 
(1997) have carried out numerical simulations of the inter- 
stellar medium in NGC 157 under the action of a bar-driven 
spiral density wave, and predicted the locations of several 
resonances accordingly. 

In an effort to see how the kinematics may be influ- 
enced by the action of a bar, we have sampled both the 
Ha and Hi velocity fields along the kinematic minor axis 
{9 = 310°), and binned the data into 15 arsec intervals, to 
match the resolution of the Hi data. As Figure ^ shows, 
both the neutral and ionised gas components exhibit almost 
sinusoidal variations about the systemic velocity along the 
minor axis. Assuming the spiral arms to be trailing, then the 
observed rotation pattern implies that the southeast side of 
the disk is closest to us. Consequently, the negative velocity 
residuals in Figure ^ along the southeast half of the minor 
axis indicate gas flowing towards the observer, and therefore 
an outflow, rather than the gas inflow more commonly as- 
sumed to be associated with the action of a bar. Thus the 
nuclear kinematics appear to be dominated by the presence 
of a mild starburst, rather than any associated bar. 
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4.3 The turnover in the rotation curve of 
NGC 157 

The extent and abruptness of the dechne in rotation veloc- 
ity in NGC 157 would, if proved, make it one of the most 
unusual rotation curves yet seen. The prevalence of flat rota- 
tion curves as seen in the ionised gas (Rubin, Ford & Thon- 
nard 1978) and in the extended Hi (Bosma 1978) has been 
taken as one of the strongest arguments in favour of the 
existence of a dark matter halo around all spiral galaxies. 
To date, only a handful of galaxies are observed to have 
truly declining Hi rotation curves, consistent with having 
reached the edge of the total mass distribution; these include 
NGC 2683 and NGC 3521 (Casertano & van Gorkom 1991), 
NGC 7793 (Carignan & Puche 1990), NGC 1365 (Jorsater 
& van Moorsel 1995), and NGC 4244 (Oiling 1996). How- 
ever, NGC 157 is somewhat more unique, in the sense that 
the decline in rotation velocity is both rapid and severe, and 
yet the outermost part of the rotation curve is still relatively 
flat, indicating that the mass distribution continues well be- 
yond the last measured point. 

There are a few galaxies where optical observations have 
also revealed a decline in the line-of-sight velocities; for ex- 
ample NGC 4303, NGC 4321 (Distefano et al. 1990), NGC 
4254, and NGC 4536 (Sperandio et al. 1995; Afanasiev et al. 
1992). However, none of these cases show a decline as abrupt 
as we have found in NGC 157. More importantly, in all cases 
cited, the decline is seen only on one side of the nucleus, and 
often the decline is not visible in the Hi rotation curve (e.g., 
NGC 4254; Distefano et al. 1990). In contrast, our Hi data 
are in good agreement with ZK's measurements, and while 
the optical decline is seen on only one side of the galaxy, 
the decline is seen at the same radius on both sides of the 
galaxy in Hi. 

Before attempting to model the mass distribution, we 
must be confident that Figure ^ is a fair and accurate rep- 
resentation of the disk rotation and orientation parameters 
in NGC 157. It could be argued that a more likely scenario 
is that the rotation curve stays flat, and the change in pro- 
jected rotation velocity seen in the outer disk is entirely 
due to variations in 6 and i. In Figure |l^, we have fixed 
Viot at 200 km s~^ (close to the maximum value attained in 
the inner disk) for r > 100 arcsec, and used the ROTCUR 
routine to again find the combination of 9 and i that min- 
imises the circular velocity residuals. These residuals average 
9.6 km s~^, only slightly larger than the 7.9 km s~^ residuals 
attained with the free fit shown in Figure ^. 

The warp behaviour of the line-of-nodes is almost un- 
changed, but the inclination is forced to drop as low as 23° 
in order to sustain such a high rotation velocity. Requir- 
ing a constant rotation velocity of only 180 or 160 km s~^ 
results in a minimum inclination of 25° and 30° respec- 
tively. The field galaxy calibration of the TuUy-Fisher rela- 
tion (Pierce & TuUy 1992) for Ms = -20.82 ±0.06 implies a 
still higher maximum rotational velocity of 234lgg^ km s~^. 
To match this velocity with a flat rotation curve would re- 
quire i ~ 35° in the inner disk, dropping to as low as i ~ 25° 
at r = 400 arcsec. Thus, it appears that the outer H I velocity 
is unusually low, rather than the inner H I velocity peak be- 
ing anomalously high, for a galaxy of this magnitude. These 
low inclinations are consistent with the axis ratios of the 
outer Hi contours (Figure H), but countering this is the fact 



that the /-band surface photometry does not support such 
low inclinations in the region 90 < r < 180 arcsec. At these 
radii, the spiral arms do not strongly disturb the isophotes, 
and to account for the difference by intrinsic ellipticity of the 
disk requires e ~ 0.2, an extreme value (c/. 0.05; Rix & Zarit- 
sky 1995). Thus, unless t he ga s and the stars are somehow 
'decoupled' (e.g.. Section 4.4.2), and orbit in totally differ- 



ent planes, the rotation curve must decline significantly in 
this region. 

One other possibility is that gas orbits in the outer disk 
of NGC 157 are intrinsically elUpttcal, violating one of the 
principal assumptions of the ROTCUR analysis, and casting 
doubt on the derived inclinations. Non-circular gas motions 
are often seen in the central regions of galaxies, usually aris- 
ing through the influence of a bar. These effects are usually 
restricted to the region near the bar (though there are ex- 
ceptions in cases where a resonance is set up between the 
orbital epicyclic frequency and the bar pattern speed: Ry- 
der et al. 1996); outside of this, streaming motions along 
the spiral arms dominate perturbations from uniform cir- 
cular motion. In NGC 157, the bar component proposed 
by Elmegreen et al. (1996a, b) and by Sempere & Rozas 
(1997) is a small-scale feature, extending to a radius of only 
~ 0.25r25 (30 arcsec), so it can be discounted as the source 



of the velocity peak in the rotation curve (Section 4.2). 

To check the effects of elliptical orbits at large radii, we 
computed a few model velocity fields and analysed them 
with ROTCUR. The elliptical orbits were not computed 
from an assumed potential, but as Lissajous figures, the sum 
of two simple harmonic oscillators. We compared orbit el- 
lipticities {e — {a ~ b)/a) of e = and e = 0.5, for a fiat 
(v) = 200 km s~^ rotation curve and a declining curve, com- 
puted from the Brandt (1960) formula for -Rmax = 80 arcsec, 
Vmax ~ 200 km s~^and n = 10. The calculations were done 
for disks inclined at i = 45° and 60° to the line of sight. 
For this small grid of models, the shapes of the outer-disk 
isovelocity contours remained almost unchanged by taking 
the orbits to be oval, but the amplitude of the radial veloc- 
ity variations was increased. ROTCUR correctly fitted the 
inclination and position angles, but for the elliptical-orbit 
velocity fields it returned circular velocities ~ 25% higher 
than the input value. For the i = 45° declining rotation 
curve case, ROTCUR had difficulty in fitting the inclina- 
tion, and did not give consistent results from ring to ring, 
or for different initial parameters. However, e = 0.5 is an 
extreme value, and given the lack of any effect in the other 
test cases, we think that the velocity field of NGC 157 can 
be well modeled by circular orbits. 

Finally, we note that although he had no knowledge 
of the extended Hi disk (nor its unusual kinematics) in 
NGC 157, Blackman (1979) did investigate the predictions 
of the spiral density wave theory by assuming a declining 
rotation curve outside of r25, of the form popularized by 
Brandt (1960). He found that such a rotation curve yielded 
a theoretical spiral pattern in much better agreement with 
the actual spiral arm pattern than did a fiat rotation curve 
in the outer regions. Although this approach involved many 
simplifications, it does provide yet another clue in favour of 
the declining rotation curve. 

In summary, the velocity field is probably not strongly 
affected by peculiar motions, but the difficulties in separat- 
ing the effects of inclination and rotation velocity remain. 
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and all we can really say is that the actual rotation curve 
behaviour almost certainly lies somewhere between the two 
extremes shown in Figures I and 0. We can then use these 
to place upper and lower bounds on the various mass mod- 
els. 



4.4 Mass Models of NGC 157 

In this section, we consider some scenarios that could 
account for the peculiar velocity field/rotation curve in 
NGC 157. 



4.4.1 'Missing' Dark Matter? 

We begin with the standard 'maximum-disk' model for fit- 
ting the declining and flat rotation curves, using the pro- 
cedure described in Walsh, Staveley-Smith, & Oosterloo 
(1997). The relative contributions to the overall rotation 
due to the stellar component, the gas surface density {am 
multiplied by 1.3 to account for the presence of He), and a 
pseudo-isothermal halo have been calculated using a least- 
squares fit to the observed rotation curve. A summary of the 
best fitting parameters is given in Table ^, and the fit to the 
rotation curve with the turnover is plotted in Figure |l^. The 
mass of the gaseous disk is 9.62 x 10^ M0, and the stellar 
disk has 1.92 x lO" M© (with a stellar M/Li of 1.4). The 
gas mass-to-stellar mass ratio is therefore quite normal at 
5%. 

In the case of the declining rotation curve, the dark-to- 
luminous matter ratio at the last measured point is unusu- 
ally low, at almost 1:1. This is almost identical to the value 
found in NGC 7793 (Carignan & Puche 1990), and only 
the massive (and distant) galaxy NGC 801 in the sample of 
BvW has a lower value. Based on the empirical relations of 
Kormendy (1990), we would have expected a central density 
po ~ 0.004 M0 pc~^ and a core radius rc ~ 19 kpc for a 
galaxy with the blue luminosity of NGC 157. This is quite 
a large core radius, and once again, only the most luminous 
{Mb > —21) galaxies in the BvW sample have larger core 
radii and lower central densities. 

Even for the case of a flat rotation curve, the dark- 
to-luminous mass ratio still only approaches 2 at the last 
measured point, though such a value is more in keeping with 
that measured in galaxies of similar Vmax and Hubble type 
(BvW). In summary, a standard gas + stellar disk -I- dark 
matter halo model points toward NGC 157 having a low 
dark matter content, though not unprecedentedly so. 

This raises the question of how the galaxy could have 
come to have a low dark matter content. Halo stripping in 
a cluster environment is easy enough to understand (Whit- 
more et al. 1988), but as noted in the introduction, NGC 157 
is rather isolated. Casertano & van Gorkom (1991) found 
that bright (massive), compact (i.e. short disk scale length) 
galaxies are more likely to show declining rotation curves 
(usually by no more than 30% of their peak velocity), though 
BvW disputes the scale length dependence. NGC 157 pos- 
sesses both a fast rotating and relatively compact disk, with 
an /-band scale length of 3.4 kpc; according to Casertano & 
van Gorkom (1991), such galaxies will have a steeper poten- 
tial well, into which dark matter is less likely to fall. Thus, 
if the dark matter content of NGC 157 is as low as the de- 



cline in the rotation curve suggests, then it presumably has 
always been low, reflecting the circumstances of its birth. 

4.4-2 Truncated or Decoupled Disks 

Casertano (1983) showed how a truncated exponential disk 
could give rise to an abrupt decline in the rotation curve just 
beyond the truncation radius. In particular, he was trying 
to account for a sudden drop in the rotation curve of the 
edge-on spiral galaxy NGC 5907, but photometry of such 
galaxies (e.g., van der Kruit & Searle 1981) had already 
made a strong case for truncations in the light distribution. 

We have used a modified version of the disk-fitting pro- 
cedure from the previous section in order to test whether a 
truncated disk could account for some (or all) of the drop in 
rotation velocity. Based partly on the findings of Blackman 
(1979), confirmed by our Figure [lo| that NGC 157 appears 
to have a second, outer exponential disk, we include 2 sepa- 
rate disk components in our models; one truncated at some 
radius Rt, and a second, less massive untruncated disk to ac- 
count for the velocity of rotation at large radii. We therefore 
found it necessary to model the bulge component separately, 
using a King profile. The scale length of the inner disk was 
taken to be 3.4 kpc, derived from a linear fit to the /-band 
profile of Figure hd in the range 35 < r < 100 arcsec. This 
disk was truncated at a radius Rt, with the best fit obtained 
for Rt = 6 — 6.5 kpc, i.e., ~ 2 disk scale lengths. To 'soften' 
the sharp fall of the density at the edge, a non-zero scale 
length (1 kpc) was introduced for R > Rt. The parameters 
of the untruncated exponential disk, which dominates the 
density a.t R > Rt were then varied to give the best fit to 
the observed rotation curve. The best fit value of the scale 
length was found to be 4.7 kpc, which is slightly larger than 
was found from the /-band surface photometry (Section |3.5| ). 
The asymptotic velocity of rotation of the model halo is only 
110 km s^^, which is about half of the maximum velocity of 
rotation of the disk. 

The rotation due to each component, plus the total com- 
pared with the observed declining rotation curve, is plotted 
in Figure |l^. The masses of individual components within 
three important radii are given in Table ^: r = 12.5 kpc, 
which corresponds to the optical radius r25 of this galaxy; 
r = 25 kpc, where the halo begins to dominate; and r — 
50 kpc, just outside the last measured point. The dark-to- 
luminous mass ratio at these points grows, from 0.13, to 0.51, 
and finally 1.56 respectively. Thus, although the truncated 
disk model can successfully account for the rapid drop in 
rotation velocity beyond 5 kpc, the outer exponential disk 
and gas components together are not sufficient to sustain 
the rotation curve, and a dark matter halo (albeit a mini- 
mal one) is still required. Taking into account the total blue 
luminosity of the galaxy (Section |3.5| ), we obtain an unusu- 
ally low M/ Lb ~ 0.6 (solar units) within r25, which further 
supports the conclusion that there is no significant dark halo 
within the luminous disk. 

Truncated disks are seldom evident in radial surface 
brightness profiles of low-inclination galaxies, most likely be- 
cause of azimuthal averaging in the analysis smearing out 
any asymmetric edges in the younger stellar populations 
(van der Kruit 1988). Radial cuts across the /3-band im- 
age of NGC 157 shown in Figure ^ give the impression of a 
more abrupt drop in luminosity towards the northwest edge 
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Table 3. Mass model properties for NGC 157. 

Declining R.C. Flat R.C. 

Parameter Vmax = 170 km s^^ Knax = 200 km s^^ 

M/Li (Mq/Lq) 1.4 0.84 0.94 

po (Mq pc-3) 0.0019 0.0071 0.014 

rc (kpc) 15 9.5 7.6 

A^dark (10" Mo)" 2.1 3.4 4.3 

"Mass of dark matter contained in halo within r < 44 kpc. 

Table 4. Truncated disk model properties within given radius for NGC 157. 

Model component Mass within R (10^ Mq) 

R = 12.5 kpc K = 25 kpc R = 50 kpc 

Truncated disk 36.0 36.0 36.0 

Outer disk 10.0 13.0 13.4 

Bulge < 1 < 1 < 1 

Gas (Hi + He) 4.5 7.1 8.6 

Halo 6.7 29.0 91.0 

Total 57.0 85.0 149.0 



of the disk than elsewhere, though this is at least partly a 
matter of contrast with the outer edge of the northern spiral 
arm. 

Truncated disks and warps might be related, since 
Sparke (1984) showed that a sharply truncated edge to the 
disk could be the crucial factor in sustaining vertical os- 
cillation modes (i.e., a warp) in the presence of a halo. Of 
course, the disk referred to in this case is the gaseous disk, 
whereas it is the stellar disk that is truncated in our model; 
furthermore, Sparke's model produces warps interior to the 
truncation. Nevertheless, the forces that give rise to one phe- 
nomenon (a warp or a truncated disk) could equally well be 
driving the other. 

We wish to stress that the phenomenon we are inter- 
ested in here is that of a truncated inner disk, rather than 
an actual edge to the total mass distribution (as might be 
produced say, by ram-pressure stripping: Gunn & Gott 1972; 
Kritsuk 1983) or to just the H I disk (e.g., photoionisation by 
the extragalactic UV radiation field: Maloney 1993; Dove & 
ShuU 1994). All these indications of a truncated inner disk 
lead us to reconsider our other data under the hypothesis 
that NGC 157 is a loosely-coupled (or perhaps even merged) 
pair of disks. The suddenness of the decline in projected ro- 
tational velocity, the coincidence of this with the optical 
edge of the disk, and the apparent discrepancy between op- 
tical and kinematic inclinations (if we require a flat rotation 
curve), all lend some support to this idea. 

The total Hi spectrum gives conflicting evidence. The 
(slightly) different systemic velocities of the inner and outer 
Hi disks would favour this scenario, but the roughly equal 
gas masses in the two components is a contra-indication, 
as the merging of equal-mass disks is expected to result in 
an elliptical galaxy (e.g., Mihos & Hernquist 1996, and ref- 
erences therein). Indeed, even minor accretions of mass are 
expected to cause significant heating of galaxy disks (Toth & 
Ostriker 1992; Quinn, Hernquist, & FuUagar 1993; Walker, 
Mihos, & Hernquist 1996). On the other hand, the Hi gas 
in the merging system NGC 520 is still recognisably 'disky', 
with a rather smooth velocity field, despite showing tidal 
tails, strong distortion of the optical light, and other signs 



of a disk-disk merger (Hibbard 1995). Similarly, NGC 4826 
is thought to have survived as a disk system after a major 
addition of mass (Braun et al. 1994). This galaxy has an 
outer gas disk that counter-rotates with respect to the inner 
gas and stars, the two disks having identical kinematic ori- 
entation parameters. Braun et al. propose that it has formed 
either by a slow accretion of retrograde gas, or the merger 
of two disks of opposite spin, one much richer in gas than 
the other. While the chances of such an orderly merger must 
be extremely slim, other examples of counter-rotating disk 
components are appearing (e.g., Jore et al. 1996), so the 
possibility that NGC 157 is also a multiple system cannot 
readily be dismissed. 



5 CONCLUSIONS 

The new radio and optical observations of NGC 157 pre- 
sented in this paper put this object among a small number 
of spiral galaxies whose photometric and dynamical proper- 
ties have been studied in detail over a large range of radial 
distances - from less than one kpc from the nucleus up to 
nearly 50 kpc, or about 14 optical scale lengths. Although 
many properties of this galaxy are rather common (grand 
design spiral structure, the total mass of luminous compo- 
nents, the maximum velocity of rotation, etc.), in some other 
respects, NGC 157 is a most peculiar system: 

(i) NGC 157 possesses a highly extended H I disk which 
is mildly warped and stretches far beyond the optical extent 
of the galaxy. 

(ii) The Hi in the inner galaxy is not only concentrated 
in the spiral arms well defined by Hll regions, but also 
shows remarkable honeycomb-like structure in the interarm 
regions, which in turn gives a strong impression of 'super- 
bubbles' in the disk. The azimuthally-averaged surface den- 
sity of H I drops sharply in the outer part of the optical disk: 
it decreases by more than a factor of three in the narrow 
range 80 — 120 arcsec (8 — 12 kpc). 

(iii) The combined velocity field of both ionised and neu- 
tral hydrogen implies an unusual shape for the rotation 
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curve, which begins to decline steeply in the outer part of the 
optical disk and continues to drop until the velocity of rota- 
tion reaches almost half of its maximum value. This points 
to the presence of a dark matter halo, the relative mass of 
which is small when compared to the mass of the stellar 
disk; a low M /Lb ratio (~ 2) supports this conclusion. It is 
argued that NGC 157 is the only known galaxy for which 
a low-mass halo may be claimed with some confidence. 

We have outlined various circumstances under which 
NGC 157 may have attained these features, including an 
unusually steep potential well, the combination of truncated 
and extended exponential disks, or possibly a merger of two 
equal-mass gas disks. Although we are not yet able to dis- 
tinguish which (if any) of these scenarios is in operation in 
the case of NGC 157, we do point out that the distinctive 
'shoulders' on the global H l profile may be useful for identi- 
fying potential analogs of NGC 157, with large and abrupt 
declines in their rotation curves. 
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Figure 1. Hi channel maps of NGC 157 at ~ 10 km s^^ intervals from our natural-weighted 'supercube'. Heliocentric velocities in the 
optical convention are given for each map in units of km s~^. The flux range plotted is from 2 to 40 mjy beam~^, and is displayed 
logarithmically to enhance details in the extended Hi disk. The cross in each panel indicates the position of the dynamical center 
(coincident with the optical nucleus). 
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Figure 2. Global H I profile (solid line) for NGC 157, obtained by summing the emission within each channel map of the natural-weighted 
'supercube'. The dotted profile is the emission from just the inner disk, while the dashed profile represents the emission from the outer, 
warped disk. 
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Figure 3. Contours of Hi column density overlaid on a _B-band image from the KPNO 0.9-m telescope. The contours correspond to 
(projected) column densities of 0.5, 1.0, 2.0, 3.5, 7.0, 10.0, 15.0, and 20.0 X 10^° cm^^. 



© 1997 RAS, MNRAS 000, §-26 



The Peculiar Rotation Curve of NGC 157 15 



-08 16 — 



o 
o 
o 

CM 



g 

I- 
< 



o 

LU 
Q 




00 35 15 



00 35 00 00 34 45 00 34 30 

RIGHT ASCENSION (J2000) 



Figure 4. Isovelocity contours plotted on a grey-scale of the same Hi distribution as shown in Figure^. The contour interval is 20 km s ^ , 
and major contours are marked. The cross marks the position of the dynamical center, as determined from the kinematical analysis. 
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Figure 5. Higher resolution grey-scale image of the uniform-weighted Hi in the inner disk of NGC 157, on which contours of the Ha 
emission has been overlaid. Notice the close, but not perfect, coincidence of star formation activity with peak gas surface density. 
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Figure 6. Radial variation of tlie depojected surface density of Hi in NGC 157, assuming a fixed inclination i = 45° and position angle 
= 220°. The 25 B mag arcsec"^ radius is marked by the arrow. 



© 1997 RAS, MNRAS 000, 0-26 



18 S. Ryder, A. Zasov, V. Mclntyre, W. Walsh and O. Sil'chenko 



o 



> 



CD 



200 



150 



100 



50 





250 



200 



20 




+ 



+ 



+ 



+ 



+ 



+ 



H — — y- 



H h 



® Both sides 
+ Receding 
X Approaching 

— \ \ \ \ K 



H h 



Ha FP 




GASP I-band 



+ 



+ 



+ 




H h 



H h 



H h 



t 



X 

+ 



¥ 

+ 



X 
+ 



+ 



_l I L 



_l I I I L 







100 



200 



300 



400 



Radius (arcsec) 



Figure 7. Rotation curve and variation of ring inclination and kinematic line-of-nodes from ROTCUR modelling of the uniform (r < 
120 arcsec) and natural-weighted (r > 120 arcsec) velocity fields. Points indicate the fits to the observed H I velocity field, the solid line 
is the equivalent result from the Ha velocity field, and the dashed line traces the orientation parameters from surface photometry of the 
7-band image. For the Hi velocity field, separate analyses have been carried out on the entire disk (•), as well as for just the receding 
(+) and approaching (x) halves. Error bars in the formal fitting are much smaller than the differences between the separate halves. 
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Figure 8. Comparison of the velocity fields obtained from Ha Fabry-Perot interferometry (grey-scale) and Hi aperture syntliesis 
(contours). For clarity, only the Ha velocity ranges 1550 ± 10, 1610 ± 10, 1670 ± 10, 1730 ± 10, and 1790 ± 10 km s"! are shown. The 
contours correspond to Hi velocities of 1550 (northeast), 1610, 1670, 1730, and 1790 km s"'^ (southwest). The synthesised beamsize of 
the Hi data is shown by the ellipse in the lower-left corner; the spatial resolution of the Ha data is 5 times better. 
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Figure 9. Contours of the uniform-weighted 1.4 GHz continuum overlaid on tlie pure Ha image of NGC 157. The Ha surface brightness 
is displayed logarithmically to bring out the diffuse component, while the contours correspond to flux densities of 1.5, 2.0, 3.0, 4.0, 5.0, 
6.0, and 7.0 mjy beam~^. 
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Figure 10. Radial surface brightness profiles for NGC 157 from ellipse fitting to the B, I, and Ha images. The disk orientation parameters 
found from the /-band fitting are shown separately in Figure M. 
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Figure 11. A cut through the uniform- weighted Hi cube at a position angle of 10°, for comparison with Figure 3 of Zasov & Kyazumov 
(1981). 
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Figure 12. Cut along the minor axes of tlie Ha and Hi velocity fields, with data binned into 15 arcsec intervals. The solid points and 
fat error bars are the mean and r.m.s. velocities from the uniform-weighted H« data, while the open triangles and thin error bars are 
from the Ha data. The dashed line is the systemic velocity derived from fitting to the H I velocity field, while the dash-dotted line is that 
from the Ha velocity field. 
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Figure 13. As for Figure ^, except that the rotation velocity is held fixed at 200 km s~^ for r > 100 arcsec. 
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Figure 14. Plot of the relative contributions to the declining rotation curve (squares) of NGC 157 from (dotted line) atomic gas, (dashed 
line) stars, and (dot-dash line) a dark matter halo, with the sum of these three shown as the solid line. 
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Figure 15. Plot of the relative contributions to the declining rotation curve (solid symbols) of NGC 157 from an exponential disk 
truncated at about 2 scale lengths, an outer exponential disk, atomic gas, stars, and a dark matter halo, as well as their sum. 
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